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ri ■ Abstract 

f^- The analytical result for the six-photon helicity amplitudes in scalar QED is presented. To compute the loop, 

O ,' a recently developed method based on multiple cuts is used. The amplitudes for QED and QED are also 

li I derived using the supersymmetric decomposition linking the three theories. 
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1 Introduction and Notation 

The light-by-light scattering is a prediction of quantum electrodynamic despite the fact that it has never been 
observed so far. The four-photon amplitudes, in QED, have been computed in the fifties at one loop for massive 
fermion [T] and recently at two loops for massless fermions in QED [2J and in A^ = 1 supersymmetric QED [3j. 
The first result for the six-photon amplitudes, at one loop in QED and for massless fermions, was obtained for 
the MHV (Maximal Helicity Violating) amplitude by Mahlon T. The complete helicity amplitudes, in QED, has 
been computed numerically by direct integration of the Feynman diagrams [51, and also, by using reduction at the 
integrand level [6 . In the same time, it has also been computed analytically using unitary cut methods and cross 
checked with a reduction method fTj. A compact formula has been given, proving the power of the unitary cut 
methods. The unitarity-cut methods were developed first in [8j and then in ^ and are currently under intense 
developments [TUl HH [HI HI E] • 

Although the six-photon amplitudes are out of reach for nowadays experiments, they provide a good laboratory 
reaction to settle efficient methods to compute one loop multi-leg amplitudes. Indeed, multi-particle processes 
involving Quantum ChromoDynamics (QCD) will play an important role in the physics probed by the hadronic 
colliders at the TeV scale. In particular at the future Large Hadron Collider (LHC), the production of four, and 
even five or six jets will not be marginal. Besides providing a refined probe of the dynamics of colour, such 
QCD processes constitute a background to the search for new particles. Indeed, the search for many of these new 
particles at hadronic colliders often relies on signatures based on cascade decays. The latter end up with final states 
involving a large jet multiplicity. Furthermore, the lowest order estimates for such processes are plagued by the 
well-known deficiencies of large renormalization and factorization scale dependencies, poor multi-jet modelling and 
large sensitivity to kinematic cuts. Therefore the calculation of next-to-leading-order (NLO) corrections to such 
cross sections is a necessary step forward. 

In this article, for scalar QED (respectively QED, M — 1 supersymmetric QED) a generic iV-photon helicity 
amplitude is denoted by: ^4.^"''"' (respectively A-'^^"^^""', A'^^'^). We will use unitarity-cut methods to compute the 
six-photon amplitudes in scalar QED. From these results, we can derive results for QED and A/" = 1 supersymmetric 
QED. To achieve this, we use a relation which relates the three theories. To find it, we proceed as follows. Starting 
with the A''-photon QED amplitudes and using the fact that degrees of freedom for internal lines can be added and 
subtracted |15j , we can write the following relation: 

^fermion r. Ascalar , /iA^=l 



A-A 



N 



2^sca(ar ^ ^VV=i (j) 



To calculate those amplitudes, we use the spinor helicity formalism developed in [1^]. For the spinorial product, 
we introduce the following notation: 

{Pa - \pb+) = {ab) (2) 

{Pa + \Pb-) - [ab] (3) 

{Pa - I Abe-) = (abc) - [cba] = (p, + | A|Pa+) (4) 

{Pa + \A i>c\Pd~-) = [aM = -\dcba\ = -{pd + Mc A\Pa-) (5) 

Moreover we use Pi...j — Pi + ... + Pj and s^. ^ — {pi + ... + Pj)'^ — pf ,-. 

The outline of the paper is as follows. In section 2, we analyse the structure of the amplitudes, and compute the 
different tree amplitudes necessary for our calculation. In section 3, we give an analytical result for the six-photon 
amplitudes in scalar QED and in section 4, we derive analytical results for QED and J\f — I supersymmetric QED. 
In section 5, we plot the different amplitudes for some kinematics and discuss potential problems. 



2 Structure of the amplitudes 



2.1 Decomposition of the amplitudes ^j*^^™""^ j^^caiar ^j^^j Aq'~^ 

The standard reduction methods, for example [T71[THl[Tn] show that any ampUtude can be written as a combi- 
nation of master integrals. This set of master integrals is not unique. 

From now, we consider theory with massless particles. In the case of the six-photon amplitudes, we use the 
following decomposition: 

Ae^ Y, («' ^4 + hF3 + d, F2A + Q F2B + e,Fi + f, Fq 

tea (1,2, 3 A, 5, e) 

+5j ^3 + "i -'3 + *!' ^3 + Ji ^2 + rational terms J (6) 

where F4 (respectively F3, Fi and Fq) is the so called "finite part" of the four point function, in n dimensions, 
with 4 external masses (respectively three external masses, one external mass and zero external mass), F2A (re- 
spectively F2b) the so called "finite" part of the n dimensional four point function with two adjacent external 
masses (respectively with two opposite external masses). Only this set of functions {-^2^,^23,^1} will be used, 
their exact definition can be found, for example in [501, and to be self consistent we recall them in the appendix [XI 
In addition, /J^mass (respectively /»2mass^ jnimass^ -g ^j^g ^ dimensional three point function with three external 
masses (respectively two external masses, one external mass) and /^ is the n dimensional two point function. The 
IR divergences are carried by the function /^^™ and /g ^™ and the UV one by the function /^. 

Using unitary-cut methods, we only have to compute the coefficients ai...ji and rational terms. Most of them are 
related by Bose symmetry or parity. 



2.2 Tree amplitudes 

In the framework of unitary-cut methods, we need to compute first tree amplitudes. In this subsection, we will 
present only the tree amplitudes in QED and scalar QED useful for our six-photon amplitude computation. The 
needed tree amplitudes, are the amplitudes corresponding to the reactions: two scalars (fermions) into A^-photons 
with the same helicity and two scalars (fermions) into A^-photons with the same helicity but one. We assume that 
all the photons are ingoing, and Pa and pb are the four momentum of the scalars (fermions): 





r igUre 1: Tree amplitudes needed for the six-photon amplitudes. The particles associated with a plain line are scalars or fermions. 



A-^(l+,...,iV+) = (7) 

4;:e™°"(i+,-,iv+) = o (8) 

N + 1 



Ascalarn+ ^r+ u~^ _ . I .. /^\' ' ' V^ ^ {kal){kbN) 

tree 



/ \ i\ i-± 

^(l+,...,7V+,A:-) = z(eV2) ^ 



.({i..A.n.)<12)(23)...(iV-liV)(lal)(iV6iV) 



.ferniion,.+ j^+ ,-^_Ascalar(^+ Ar+ U- \ f V^ , itzL] 

^tree l-L 7 •••7 J* 7'* ) ~ ^tree l-L 7 ■■•7 J* 7 *« '\(Uh\ (hfj\j 



(9) 
(10) 



where Si — —eV2- — -— — is the eikonal factor. Note that in equation pU|) . a sum over the heUcities of the fermions 

{ai){ib) 

has been performed. 



2.3 Additional reductions 

Now using the properties of the QED theories, we can simphfy furthermore the decomposition ^ of the six- 
photon amphtudes. Those comments and rules is available only for these amplitudes. 

Remark 2.1 Only the functions I^ '",^3 ™ are IR divergent. Since each diagram is not IR divergent, the coeffi- 
cients hi and ii are zero. 

Each Feynman diagram of the six-photon amplitudes is free of IR divergences thanks to the numerator of the 
fermionic propagator in QED or the structure of the vertex in scalar QED. If the reduction is done by pinching 
propagators, we get sub-diagrams which are not IR divergent. After the reduction, we obtain three point sub- 
diagrams and the "finite" part of four point scalar integrals. Since the three point sub-diagrams are free of IR 
divergences, they cannot be expressed in term of one mass/two mass three point scalar integrals I^'^"^ , /g ^"^ and 
so the coefficients hi and ii are zero [T5] . 

Remark 2.2 Using standard reduction (for example 119]), we can show for QED and scalar QED, that the coeffi- 
cients in front of two point functions are zero and also that the rational terms are zero. 

The first statement is in accordance with the fact that each diagram of the six-photon amplitudes is free of UV 
divergences. From that, it is not obvious that the different coefficients are zero, we have proven it by explicit calcu- 
lation. As a consequence of remarks (12. ip and (|2.2p . there are no logarithmic terms in the six-photon amplitudes. 
The second statement, shown in ref. [21| . is also far from being obvious. Indeed, according to power counting 
arguments |22j . the rational terms can be present. In fact, they are present for individual Feynman diagram but 
these rational terms sum up to zero when adding all the diagrams. 

Now, from what has been said previously, we can reduce the decomposition of the amplitude © : for the six- 
photon case, the coefficients hi,ii, and ji are zero and the rational term is also nul. In addition, since we have only 
six photons on shell, the coefhcients ai, hi and fi are also zero. So each amplitude ^J'='''"*°"^ ^^caiar ^^^^ ^g written 
as: 



ieCT(l,2,3,4,5,6) 



d^F2A + c,F2B+e,F^+gaT""'' (11) 



From the helicity structure of trees, we can derive some rules which will reduce furthermore the decomposition ()lip . 

Rule 2.1 Consider a master integral with mass. If the mass is formed only with photons with the same helicity, 
therefore the coefficient in front of the master integral is zero. 

Proof : Using the cut technics, we get that the coefficient in front of this integral is proportional to the tree 
amplitude which, once pinched, yields the mass. Formulae (I7I8|) show that the on-shell tree amplitudes with 
photons having the same helicity are zero. So the coefficient of a master integral with a mass formed by photons 
with the same helicity is zero. 

Rule 2.2 For the box with one mass and with two adjacent masses, the helicity of two adjacent massless legs must 
be alternate. In the case of the box with two opposite masses, the heUcities of the two opposite massless legs must 
be the same. If it is not the case, the coefficient in front of the master integral is zero. 

Proof : In the case of the one mass box and the two adjacent mass box, if the helicities of two adjacent massless 
legs are the same, the coefficient, in front of the box, will be proportional to trees, eqs. (|7l8p . which are zero. The 
case of the two opposite mass box is more complicated and a proof is given in the appendix [BJ 

Now we can further reduce the decomposition of all helicity amplitudes thanks to these last rules. We begin 
with the most simple amplitudes: AAr(l^, 2+..., Af+). With at most one negative photon, we can have only one 



mass according to the rule (|2.ip . so Ci,di,gi ~ 0. But in this case, the hehcities cannot be alternate so, from the 
rule (|2.2p . we deduce that each e^ = 0. Therefore, we get for these amplitudes: 

VA^ > 4, ^^""''^'^(1=^, 2+..., N+) = (12) 

VA^ > 4, ylJ^"™""(l±, 2+..., N+) = (13) 

VA^>4, A^=i(l±,2+...,7V+) =0 (14) 

These results were already found by Mahlon [4 many years ago. 

Then we study the MHV amplitude ^g(l~, 2~, 3+, 4+, 5+, 6+). We have only two negative photons so we can 
only have at most two masses: gi — 0. The hehcities must be alternate according to rule (|2.2p therefore for the 
MHV amplitude, di = 0. So we have: 

^ei 1- + + + ) = ^ , g 2 X! X! -^-P'2b(s215,S415,Si5,S26) + -^Fi(s23,524,Sl56) (15) 

cr(l,2) cr(3,4,5,6) 

Some permutations leave master integrals invariants, that is why we divide the coefficient by the adequate number. 

It is rather easy to show that, in the case of the MHV six-photon amplitude, the two coefficients di and e^ are 
related. Indeed, we can take the limit such that the photon with hclicity "+" forming one of the mass of the 
two opposite mass box is soft. In this limit, wc have the following relation: lim„i2^o^2B = ^i- Since the MHV 
amplitude is composed by two MHV trees already expressed in term of eikonal factors eq. ^ and since the five 
photon amplitudes are zero, we can deduce that Ci = — e^. Finally, we have the following decomposition for the 
MHV amphtude: 

Aei 1- + + + ) = » , p 2 X! X! -^ (^2b(s315,S415,Si5,S26) -^l(s23,S24,Si56)) (16) 

cr(l,2) cr(3,4,5,6) 

We have to compute only one coefficient. 

Lastly, we examine the form of the decomposition of the Next to MHV amplitude (NMHV) Aq{1~ , 2~ , 3~ , 4:~^ , 5~^ , 6~^) . 
Here we can have three mass three point functions. But the rule (|2.2|) imposes Ci = 0. So we have only three coef- 
ficients to calculate and the amplitude can be written as: 



^ei ^+ + ) = i ^^^^ ^ ^ 2 dii^2A(si4,S452,S25,S36) + -;^^l(s63,S61,S425) 

(t(1,2,3)(t(4,5,6) 



e Q 

+ -TrFi{s25,S24,Si36) + -^Is^"" {su, S25, 835) (17) 

2 D 

Here again, we divide the coefficients by the adequate number to take into account permutations leaving invariant 
the corresponding master integrals. There are two kinds of one mass box integrals. The first kind has two photons 
with a positive hclicity and one with a negative helicity forming the mass, whereas the second kind has one photon 
with a positive helicity and two photons a with negative helicity forming the mass. As we have three photons with 
a negative helicity and three with a positive helicity, the two kinds of one mass box integrals are directly related 
by parity. That is why the coefficient in front of them are complex conjugate. 

In the next section, we will give the results for the various cases. 



j^^caiar amplitudes 



3.1 ^^'^'"'■(1-, 2", 3+..., A^+), A^ > 4 helicity amplitudes 

In this section, we calculate the MHV amplitude ^^"''"■(l-,2-,3+...,7V+) for N photons with TV > 4. The 
generalization to TV photons is rather easy, but note that we have to consider A^ > 4 otherwise there are some UV 
problems. Using the quadruple cut technics [3], by direct computation, we obtain: 

(p^/o\ ^scalar 

ct(1,2)(t(3..JV) ^ ' 

+ * 167r2 E E E fTV - MV(M ~ 4V-^2s(si35...Af , S145...M, S15...M, S2M+1...n), (18) 

cr(l,2) a{3..N) M=5 ^ ^'^ ^' 

where 

.scalar (34)^-^ (13) (41) (23) (42) [34] 



The factorial coefficient {N — M)l{M — 4)! and {N — 4)! are the number of permutations which leaves invariant the 
mass of the master integral. The formula obtained is explicitly invariant by exchange of the two photons with a 
negative helicity. In the case where A^ = 6, the amplitude reduces to: 

.scalar ( , , , .^ ■ {^^^ y ,-^ (13) (41) (23) (42) m F,is23,S24,S,,e) - F2BiSl35,SU5,Sl5.S26) 

' ^ ' 167r2 Z. Z^ (35) (45) (36) (46) (34) s^ 

cr(l,2) cr(3..6) 

(20) 

It remains only one Gram determinant (S34) in the denominator but _Fi, F2B — sai when S34 — > 0. Therefore the 
potential numerical problem when the Gram determinant vanishes is under control. 

3.2 A^^'^'"^(l-,2-,3-,4+,5+,6+) helicity amplitude 

We use the quadruple-cut technics to calculate the box coefficients [9] , and the triple-cut technique to calculate 
the triangle coefficients [23]. We obtain for the amplitude: 

(e\/2] pScalar 

^s^calar^ + + + )=,1_^ ^ ^ 2 df" ""-^2^ (si4, ^452 , S25 , S36) + ^^— i^l (^63 , S61 , S425) 

fj(l,2,3)CT(4,5,6) 

pScalar * ^scalar 

+ ^^^ Fi(s25,S24,Si36) + ^^-^/3"""(si4,S25,S36) (21) 

2 D 



where 



5425(24) [16] [Ip4252][6p4254] 

(45)[31][lp4254]2[lP4255][3p4254] 

,;„, (2P4251)(2P4253)[36][16]S425 (31) 



(22) 



scalar \^i:"izo^i \-'t"i^o^l l""J L^"J"'i^o \"^/ (')'\'\ 

(4P4251)(5P4253)(5P4251)(4P4253)[31]S31 ^ ' 
^scalar* ^ [6^6135] [6^6134] (42) (52) 5513 [45] 

[Ip6134][lp6135][3p6134][3p6135](45)s54 ^ ' 

scalar ^ [4p25l] [5pi42] [6^253] y [1J^2^1] [2J^|2] [3/^2^3] 

^' [1P254] [2pi45] [3p256] ^ [^KlA] [5X|5] [&K^2^] 

Kf = 7± (-P25)^ - S25 (P14)'' (26) 

7± =-P25.Pi4±\/A (27) 

A = (P25.P14)'-P?4P25 (28) 



The result is very compact. By expanding the coefficient g^'^'^'-"-'^ ^ we find that all square roots, coming from 7±, 
disappear as it should be. All the coefficients are rational functions of spinor products of external momenta. There 
is also a Gram determinant in the denominator of the coefficients (i|^°''"" and g^'^"-'-"-^' . These two Gram determinants 
go to zero in the same phase space region. For this region, in the numerator, it is a combination of F2A and Jg 
which cancels in such way that there is no singularity (for more details see [E]). We do not give here the explicit 
formulae because they break the simplicity of the expressions but we implement them in our numerical code. It is 
particulary important to do this for the scanning of the Landau singularities. 

In the next section, we present the result for the amplitudes with a fermion loop and a sfermion loop. 



Afermion 
6 



and Aq ^ amplitudes 



4.1 ylj^"™°"(l-,2-,3+...,A^+) and A^=\l- ,2' ,3+..., N+) helicity amplitudes 

To calculate the MHV amplitudes : Aj^'=™""(l-, 2", 3+..., N+) and Aj!^=^ {1- , 2- , 3+ . . . , N+), we use extensively 
the formula (fT(]| at the integrand level. The idea is the following. The introduction of (fT(7| under the integral, 
in the cut fermion amplitude Ay^^^^°" ^_^_^, allows to write the cut fermion amplitude as the cut scalar amplitude 
j^scaiar ^^^ plus somc terms. Using the supersymmetric decomposition ([1]) the remaining part is identified with the 
cut supersymmetric amplitude A-j^^^^^^. So, we do not have to calculate all the supersymmetric diagrams. 

We show how it works on an example. We consider a fermion loop with two cuts on propagators qa and qt- The 
figure [2] shows how the helicities are shared: two trees with one photon with a negative helicity. We compute the 



chain 2 




chain 1 



Figure 2: Spinor qed loop. 
cut amplitude ^^"^^"^cut corresponding to figure [2j 
<"™°"cut = -JdyS (ql) S {ql) {b (n+ + n_) (chainl)a>(a (n+ + n_) (chain2)6) 



dy6{ql)S{q',) 



[lb) (2a) {la)(2b) \ 
(la) {2b) {lb) {2a) J 
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scalar 



A' 



scalar 



d^q5{ql)5{ql)A{ 



scalar 



A' 



scalar 



{12)' rqS{ql)S{q!) 



t.2 Ascalar Ascalar 

2\ X /„2\ *afc^l tree^2 tree 



{labl){2ab2) 



(29) 

where n± — —^ are the chiral projectors. We recognize the cut scalar amplitude A|^°''"'^^j in the first term 

of the right hand side of equation (|29p . Since the cut amplitude A^^™''^°^\^i also obeys to the supersymmetric 
decomposition (jT]), that means that the second term of the left hand side of eq. ([^5]) is identified as the cut 
supersymmetric amplitude. Using this trick, we can easily calculate the supersymmetric amplitude and obtain 
straightforwardly the spinor amplitude. For these two amplitudes, we get: 



A 



f ermion / J\f^l 
N 



(eV2) 



N 



E E 



jfermion/Af—l 



167r2 '^ '^ (TV -4)! 

cr(l,2) cr(3..JV) ^ ' 



--Fl(s23,S24,Sl5...Af) 



+ i- 



3^2) 



N 



167r2 



cr(l,2) ct(3..A')M=5 



W-1 / ,NM-6j/e''™«o"/^=l 

(7V-M)!(Af -4)! 



-F2B(si35...Af , S145...M, S5...M7 SM+l...Af), (30) 



,fer^^on _ , (34)^-^ (13)^(42)^ 34] ^^^ ^j^^, ^ (34)^'^ 

nf=5(3j'>(4j> ^34 (34)' ' nf=5(3j>(4j) 



wnere d^ = 2 — j^^ ToTT' ^'^'^ "'i = ]v — ^ (12) . In the case or six photons, we 



get: 

.ferrmon, . (e\/2) ^ ^ 2(13)^(42)^ [34] f 1 (S23 , ^24, 5l56) " ^25 (S135, S145, Sl5, ^26) 

^ ^ ' 167r2 Z. Z^ (35) (45) (36) (46) (34) S34 

cr(l,2) cr(3..D) 



(31) 

A^=\- - + + + + ) = I , X] II 7^7w^^w5gw^(J^l(s23,S24,Sl56)-i^2s(si35,Si45,Sl5,S26)) 

cr(l,2)cr(3..6) ^ '^ '^ '^ ' 

(32) 

Full agreement is found with [HIT] for the fermion amplitude ^^<='''"'°"( h + ++). We can do the same comments 

as in scalar QED. We point out that the Gram determinants S34 have disappeared in the supersymmetric amplitude 

A^=H- - + + ++) • 

4.2 yl{'^"™°"(l-,2-,3-,4+,5+,6+) and A^=i(l-, 2", 3-,4+, 5+, 6+) helicity amplitudes 

To calculate those two amplitudes, we use again quadruple cut technics [9] and triple cut technics [23]. As in 
the preceding subsection, we use the supersymmetric decomposition ^ to extract the supersymmetric amplitude 
and then obtain the fermion one. 



Figure 3: Spinor qed loop. 

As in the previous subsection, we treat an example to illustrate how it works. It is a bit different here because the 
trees do not have the same helicity structure. We consider a fermion loop with two propagators cut as shown on 
figure [31 If we compute the corresponding cut amplitude, we get: 

<'''"""c„t = - j<r\5 {ql) 5 {ql) {b (n+ + n_) (chainl)a)(a (n+ + n_) (chain2)6) 

//• A scalar A scalar 

d-qS (ql) S (ql) Ar'^\^^Ar"'\ree " (1^2636)^ J d^qS (ql) d (ql) ' ^/g^/^g^ "^^ 

(33) 

Again, thanks to the supersymmetric decomposition ([l]), we can identify the supersymmetric cut amplitude and 
obtain the fermion cut amplitude. Doing all the calculation, we get for these two amplitudes: 

ct(1,2,3) cr(4,5,6) 

C/erTnzon/A/'— 1 
2rf- -f2A(Sl4, 5452, 525, 536) H ^ ^1 (^63, ^61, 5425) 

fermion/J^—1 * fermion/ J\f^l \ 

+ ^^ ^ i^l(s25, S24, S136) + ^^ ^ /3"'"(S14, ^25,536) (34) 



where: 



.ferrmon ^ _ 1 [1^4252]^ [6^4254]^ + g^ (24)' [16]^ 

[31](45)[lp4254]2 [1P4255][3P4254] ^ ' 

fermion __ „ (2^4258) [16] S425 (31) , . 

^' ~ (4p425l)(5p4253)(5p425l)(4p4253) [31]s31 ^ ' 

fermion* _ ^ [6^6135]' (42)^3613 [54] 

[1P6134][1P6135][3P6134][3P6135](54)S54 ^ ' 

gjerrmon ^ _^^scalar (gg) 



fermion _ ty seal 



Again, we get a full agreement with [7J for the amplitude ^j^'''"*°"( 1-++). We can observe a factor "2", except 

for the coefficient of the two adjacent box between the coefficients of the two amplitudes yiJ'='''"*°"( 1_ ^_|_) 

and Aq'^°'''"-^{ h ++). This factor "2" comes from the fact that for a fermion loop, there are two currents. 

These two currents give rise to a factor "2" except in the case of the two adjacent mass box, where they give rise 
to a sum two terms (eq. ([35])). 

dA^=i ^ I [^P4252]' , , 

[31] (45) [Ip4255][3p4254] ^ ^ 



e 



X=i _ o^.V=i 



e; 



= 2df =1 (40) 

f-'* = ef-' (41) 

gf-' = (42) 

In the Af — 1 supersymmetric QED, the coefficients are simpler. Since there are only six photons, the two kinds of 
one mass four point boxes are related by parity. This leads to the equality between the coefficients e^^^ — ef"^^. 
With eight photons or more, this becomes wrong. We can also point out that all Gram determinants disappear. 

4.3 Absence of triangles in A^^^ 

In the JV = 1 supersymmetric amplitude Aq^^, the fact that there are no triangles is probably just an accident. 
But we can make the following statement: for one loop, in QED-^^^ theory, the N photons NMHV helicity 
amplitudes will have no triangle function. In fact, as we have only three photons with a negative helicity, each mass 
of the triangle is formed with one photon with a negative helicity. The relation (|10p shows the linearity between 
the MHV tree amplitude in QED and in scalar QED. So we obtain directly using [23] that: 

H—l ferjnions , n scalar (A'^\ 

But now if we have at least four photons with a negative helicity, one mass of the triangle will be formed by two 
photons with a negative helicity and the tree amplitude corresponding to this mass will not be a MHV tree. The 
problem comes from the fact that the relation linking NMHV tree amplitude in QED and scalar QED is not linear 
but affine. This affine coefficient is the contribution of the triangle to the supersymmetric Af = 1 amplitude and 
there is no reason that the coefficients in front of the triangles become zero. So we can conjecture that at one loop, 
in QED^=^ theory, the iV > 6 photons Next to Next to MHV (NNMHV) helicity amplitudes will have triangles. 



5 Numerical results 

In this section, we present numerical results for the different six-photon amplitudes. We have built a fortran 
code using the GOLEM library for the different four point and three point functions. To be consistent with previous 
results, we use the kinematics defined by Nagy and Soper in ref. 5'. First of all, we recall this kinematics: the 
reaction 7'^^(fci) + 7'^^(fc2) -^ 7^''(fc3) + 7'^''(fc4) + 7^^(fc5) + '^''^'^(ke) is considered, where ki is the four-momentum 



of the photon i and A^ its hehcity, the four- momenta fulfiU ki + k2 = k^ + k4 + k^ + fee- In the center of mass frame 
ki + k2 = 0, fci is along the -z-axis, an arbitrary phase space point is chosen: 



k3 = (33.5,15.9,25.0) 
kl = (-12,5,15.3,0.3) 
A^ = (-10.0,-18.0,-3.3) 
h^ (-11.0,-13.2,-22.0) 



(44) 



Then new final momentum configurations is generated by rotating the final state through angle 6 about the y-axis. 
For all the plots of this section, we take a = e'^/Air = 1. The helicities "+" and "-" refer always to ingoing photons. 

5.1 The MHV amplitudes. 



In the figure m we plot the module of the MHV amplitude for QED, scalar QED and J\f = 1 supersymmetric 
QED (respectively the formula ((3T|) , ([20|) and ((32|) ) against the variable 9. Note that these formula have been 
derived assuming that all the photons are ingoing. In order to match previous results of the references [5], [7], 
we compute y^J'^""'°"/^'='* '"'/ - (/j^^ /j;2, — fcs, — fc4, — fcs, — fcg) with the helicities Ai = — , A2 = — , A3 = +, A4 = +, 
As = +, As = +. 



45000 




r igure 4: The MHV six-photon amplitudes with the Nagy and Soper |5f configuration for the three theories. Note that the curve for 
^sca ar Q^jYiplitudc has a minimum which is not zero. 

All these MHV amplitudes are tt periodic. 



5.2 The NMHV amphtudes 



In figure [3 with the same configuration than in the last section, we plot the NMHV six-photon amplitudes 
for the three theories against 9. We compute y^J^''™*°"/*^'' '^''/ - (^/j^, — ^3, — /cg, fci, — A:4, — fcs) with the helicities 
Ai = +, A2 = — , A3 = — , A4 = +, A5 = +, Ag = — . 

In this case, the amplitudes are not tt periodic. The dips appearing in the curves, are related to the Landau 
singularities called the "double parton scattering" ;5j. 



5.3 Double parton scattering 



The Landau equations give the necessary conditions for a Feynman diagram to have a singularity. In the case 
of the six-photon amplitudes, since all internal and external particles are massless, three types of singularities can 
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r igUre 0! NMHV six-photon amplitudes with the Nagy and Soper TE^ configuration for the three QED. 

appear. Two are the well known soft and coUinear singularities, the other corresponds to the so called double parton 
scattering. In this case, four propagators are one mass shell, these propagators are adjacent by pair. 

So here we explain rapidly what is the "double parton scattering" kinematics. 




Pi 



P3 P5 

r igUre D! Double parton scattering configuration: pi,p4 are ingoing photons and P2,P3,P5,P6 are outgoing photons. 



The two ingoing photons 1 and 4 split into a fermion anti-fermion coUinear pairs, each fermion scatter with an 
anti-fermion to give a photon pair with no transverse momentum in the center of mass frame pi +Pi — (c.f. figure 
[H|) . In the configuration of " double parton scattering" , the two propagators qi and q2 are coUinear to the external 
leg pi and the two other ^3 and q^ are coUinear to ^4: 



qi = -xpi 
92 = (1 - x)pi 

^ 94 = (1 - y)p4 

Solving the Landau equations, we find that the conditions to have a double parton scattering singularity are: 

det{S) -^ 

S35,S26 > 
Sl35,S435 < 

with 

det{S) = S135S435 - S35S26 



(45) 



(46) 



(47) 



In the center of mass frame pi + p4 ^ 0, det{S) — S14 fc| where k^ is the square of the transverse momentum of 
the photon photon pairs 2,6 and 3,5. 
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In the fig. [71 we plot the NMHV QED amphtudes as a function of 9. On the top of that, we surimpose k'^, 
normalized in such way that the curve is visible. We note that the dips appear in the region where k^ is minimum. 
In this case, the minimum of k^ is different from zero because we are not sitting on the singularity. A more detail 
study, will be presented in a forthcoming publication [2i] . 
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r igure 1\ Localisation of a Landau singularity. Note that the red curve describing the k^ of the photon pairs 2,6 and 3,5 does not reach 



6 Conclusion 



In this paper, we have obtained all six-photon hclicity amplitudes in QED, scalar QED and M = 1 super- 
symmetric QED. Those amplitudes are linked among themselves by the relation ([T]). To calculate them, we used 
the powerful unitarity-cut technics and we got very compact expressions. More work is required to understand 
quantitatively the behavior of these amplitudes and especially if the numerator can regularize the double parton 
scattering singularity in a kinematics where it shows up, this will be presented elsewhere |24j . 
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A Scalar integrals 



In this appendix, for sake of completeness, the definition of master integrals used in this paper is recalled, more 
details can be found in [20]. We also give det{G) the determinant of the Gram matrix Gy = 2pi.pj built with the 



external four momentum and det{S) the determinant of the kinematical S-matrix defined by Sij 
the Qi arc the four momentum flowing in the propagators. 



[I J 



where 



A.l Three mass three point function 



pi 



Pi 



P3 



Pi P5 



Si4 
S25 
S36 



11 



j: 



3 [mi,m2,m^) 



2Li2 1 - 



2 xi 



y2 



2Li2 1 - 



X2 



"3" 



where: 



xia = 



yi.2 






ml + ni\ — TO3 ± VA 



2 /2;2 



y2 



2 / 2^2 



In- ri + in^ n 1 In- :i± 



2/1 



2 / a^i 



y2 



2mf 



^l — to| + m| ± VA 
2mf 



/I /I /I OQ OQ QO /Q\ 

A = TO^ + ?7i,2 + TTig — 2?77,]^?Ti2 — 2?7i]^TO3 — 2m2m^ — i sigTi(mi) e 
The formula (jA.l[) is vaUd in all kinematical regions because of the small imaginary part i e: 

VA±ie, A> 
±i^f^, A < 



VAiie = 

The two determinants are given by the relations: 

det{G^rn) — mim2 ~ {rni.m2Y 
det{S^m) = 2m\m\m^ 

A. 2 Four point functions 
A. 2.1 With zero mass 



A 



P2 



-Vi 




P-i 



12 is,t) =^^ {(-.)- + i-t)-^} - ^Fo(.,i) 



2 

'st 



where: 



Fo(M)^i{<n=(f)+>r^} 

The determinants are given by: 

det{Go) = -2st{s + t) = 2stu 
det{So) = (st)^ = (24342)2 



(A.l) 



(A.2) 

(A.3) 
(A.4) 

(A.5) 



(A.6) 

(A.7) 



(A. 



(A.9) 



(A.IO) 
(A.ll) 



A. 2. 2 With one mass 



pi 



P:i 



, P6 

-Pr, 
' Pi 



s 


:^ 


SI2 


t 


= 


S23 


u 


■^ 


Sl3 



S456 



12 



rr 



12 {s,t,m^) = -ij {((-.)- + H)-) + ((-.)- - (-m^)-) + (H)- - (-m^)-)} 



ste 



st 



Fi{s,t,m^) 



Fi (s, t, m^) = Li2 U - ^ ) + i«2 U - ^] - ^^2 (-^) - Li2 (--^ 



where: 



= Fo{s,t) + lLi2\l-—\ +Li2\l- — ]- — 

The determinants are given by: 

det{Gi) = -2st (s + i - m^) = 2siu 
dei(S'i) = {stf = (2TO3m2)2 

A. 2. 3 With two adjacent masses 



P3 



t P5 



Pi P4 



s = 


Sl4 


i = 


S425 


m^ 


= S25 


777-2 


= S36 



(A.12) 

(A.13) 
(A.14) 



(A.15) 
(A.16) 



/" {s,t,ml,ml 



rr 



{{~s)-^ + (M)-^ - (-".?)-) + ((-i)- - (-m^)-)} 



(st)e2 



where: 



(A.17) 



/ jY) \ I jTi \ 1 / S \ / Til \ 1 

F2A (s, i, m?, 777^) = Li2 ( 1 - X ) + ^'2 ( 1 - X ) + 2^^" \V ^"^ vf) "^ 2^ 
The determinants are given by: 

det{G2A) — ~2s {ra^ni2 — t{m^ + ttij — s — t)) = — 2s(l777i477i2l) 
det{S2A) = (st)2 



s \ ( m 
n I — rr in 



A. 2. 4 With two opposite masses 



P2 t 



Pi 
P-l 



■ PC, 
>5 



P3 



S = Sl43 

t = S243 

< W = S23 

771^ = Sl4 
l^ 777^ = S56 



(A.18) 



(A.19) 
(A.20) 



12{s,t,mj,ml) 



rr 



{st — m\m'^)e^ 



{((-.)- - (-777?)-) + ((-.)- - (-7772)-)} 
{ ((-t)- - (-777?)-) + ((-i)- - (-7771)-) } 



[st — m\m'^)e^ 

— 2~^^2B (s, t, mj, ml) 

st — 777|7772 



where: 



(A.21) 
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F^B{s,t,m\,Tr?^ = -iz2 ( 1 ^j + iia ( 1-^ 

2 \ / ,„ 2 



Lz, (l - ^) + i..2 ( 1 - f ) + ^^2 ( 1 - ^) + ^^-^ (?) (A.22) 



9 2 ? 



= Fi (s,i,mf) +i^i (s,i,m^) - F^{s,t) - iLi^ f 1 - -^ - — \ \ (A.23) 

The determinants are given by: 

det(G2B) = -2 (mim2 - st) {rn\ + ml - s - t) ^ 2u {st - mlml) (A. 24) 

det{S2B) = (si - mlmiy = (2mi3mi2)2 = (2m23m22)2 (A.25) 



B Proof of rule ( 12:21 ) 
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In this appendix, we want to prove that, in the case of the box with two opposite masses, the hehcities of the 
two opposite massless legs must be the same otherwise the coefficient in front is zero. To do that, we consider the 
foUowing box integrals where the helicities of the two opposite massless legs is different. 

yi e [f ..6], pf = 

P23 ^P2+P3 
P56 ^P5+P6 

We assume that the helicity of the photon pi is positive and the helicity of the photon p^ is negative. Using the 
four cuts technics, the coefficient, called C, in front of this master integral is given by, in scalar QED: 

C (X Y^ et .que:^ ■<i2z (B.26) 

where qi^ and q2i are the solution of the four cuts conditions: 

5{ql) = (B.27) 

5{ql) = (B.28) 

5{ql) = (B.29) 

5{ql) = (B.30) 

So first we solve this system and after we will calculate (|B.26p . 

We choose as a base of the four-dimension Minkowski space: B = {^1,^4, (17^4), (47^1)}. In our case, q^ can 
be taken as a four-dimension vector, therefore, we can project it on the base B: 

q'i^ap'l + b p^l + ^(l7'^4) + ^(47'^1) (B.31) 

So to know the vector q^, we have to calculate, the four coefficient a, 6, c and d. The conditions (|B.27p and (|B.30p 
impose : 

(gi -pi)2 = o^2(pi.gi) =0^6 = (B.32) 

The conditions ((R28l) and (|R29)) impose: 

(gi +P23 +P4)' = ^ 2pi. {qi +P23) = ^ asu + 2(P23.P4) = ^ a = -^l^H^ (b.33) 

Sl4 
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The first condition (|B.27p . knowing 6 = 0, imposes: 

g2^0^cd = 0^c = 0ord = (B.34) 

If we assume that d — 0, therefore the conditions ()B.27|) and (|B.28p impose that: 

{qi + P2zf - ^ 2(p23.gi) = -S23 ^ c = ^^E^ (B.35) 

^ fl 11 Kf ■ M (4^237^14) 

and finally we obtam q^ = . 

2S14 

Else if we assume c = 0, therefore the conditions (JB.27P and (|B.28[) impose that: 

(gi + P23)' = ^ 2(p23.gi) = -S23 ^d= ^i^^ (B.36) 

Sl4 

[4^237'' 141 
g^ = . Finally according to the four cuts technics, the loop momenta is 

2si4 

^ (4^237^14) 

? ' " [4P23yi4] (B.37) 

lib = TT 

From the formula (|B.37p . we can compute q2a/b ^'^d we obtain: 

M M ^ M (47^^2314) 

92 a = lla + P23 - 



found to be: 



(l2b = 1lb+P23 



[47^'P23l4] (B-38) 

2si4 



We are now ready to compute the left hand side of equation (|B.26|) inserting the formula (IB.37|) and (jB.38|) . we 
obtain directly that: 

C ex e^.qi^e^.q2a + ei .qii^e^ ■q2b = (B.39) 

Therefore the hypothesis that the two photons pi and p^ have two different helicities implies that the coefficient in 
front of the two opposite mass integrals is zero. 
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